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Abstract: In this paper we present a general formula for the inhomogeneous non-Gaussian 
integral IdiSi,S2) = J dxi . . . dxde~2^i~^^ ^ where Si and 52 are symmetric quadratic 
forms. The solution depends on the eigenvalues of the matrix A = —iM2M^^, where Mi 
and M2 are the matrix representations of Si and 52 respectively. In the 2-dimensional 
case we also give a manifestly S'0(2)-invariant formulation in terms of invariants of the 
matrix A. An expression for I{Si, S2) in the infinite-dimensional case is calculated and the 
solution depends only on the determinants of Mi and M2. The infinite-dimensional case 
may be of use in QFT. 
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1. Introduction 

Non-Gaussian integrals with homogeneous symmetric forms of degree three and higher 
were recently studied in [P and several low-dimensional cases were solved. In this short 
note we investigate the non-Gaussian integral 

Id{Si,S2)= j dxi...dxde-"^^l-^^, (1.1) 

where Si and S2 are symmetric quadratic forms. This is the simplest case of a non- 
homogeneous, non-Gaussian integral. The problem is simplified by the fact that both Si 
and 5*2 are quadratic since, as such they can be expressed in terms of matrices. We therefore 
rewrite the integral Id in the following form 

Id{Si,S2) = j dxi...dxde-'^^'-^^ 

dxi... dx^e-i^^'^^-^i^^'-^^^^a; 



dt , (1.2) 

yJdet{M2 - it Ml) 



where we have used the well-known result 



dxi...(ixrfe-^^^^ = ^^=. (1.3) 
y/detM 
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In Q, A. Morozov and Sh. Shakirov calculate the integrals in terms of the S'L ((i)-invariants 
of the homogeneous symmetric form S. As it turns out that approach is not the most 
convenient for the present problem. Instead, we shall rewrite the integral to become a 
function of the eigenvalues of the matrix A = —iM2M^^, in the following way, 



Id{Sl,S2) = j dt 



y/det{M2 - itM_ 



at 



^det{Mi) J ^Jdet {tl - A) 



d/2 f fi-*V2 

dt- 



^det{Mi) J y/{t- zi)...{t- Zd) 



^Jdet{Mi) 

where Zi is the z:th eigenvalue of the matrix A and 



Jd, (1.4) 



„f2 12 

= I dt \ =. (1.5) 

V(t -Zi)...{t- Zd) 



2. Ward identities 



In this section we investigate differential operators that annihilate Jd- Because of the high 
symmetry of the integrand they take a rather simple form. For simplicity, let us consider 
the one-dimensional case, where 

J.=y.*^. (2.1) 

We would like to find an operator which turns the integrand into a full derivative. To that 
end, consider 

dJi 1 r -ze-*'/2 
-z—^ = - dt- 



dz 2 J {t- z)3/2 
- Ji - - / dt- 



1 3 r e-*V2 
-^1 + T / dt- 



2 4 7 (t - z)V2 

1 d'^Ji , , 



which is equivalent to 



TT + 2 

dz^ dz 2 



dJi 1 

— i + - Ji = 0. (2.3) 



We have so far not specified the contour of integration in equation 2.1. It is crucial in 



the derivation of the Ward identity that the boundary terms vanish. To that end we only 
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consider closed contours. In analogy with we denote an admissible contour as a contour 
for which the integral |2.1| converges. In this case any contour that tends asymptotically to 
the lines Arg{z) = and Arg{z) = vr are admissible. To make the integrand single-valued 
we introduce a branch cut along the line Arg(z) = as indicated in figure |^. Note that 
the contours shown in figure |l| are closed on the Riemann sphere if the infinitely remote 
point is taken into account. There are two classes of admissible contours, and these 



,nrrLii(li nit 




Figure 1: Contour of integration. 

classes are separated by the branch cut. The fact that there are two essentially different 
contours means that there will be two different solutions to the integral 2.1. This can also 
be seen from that fact that the differential equation 2.S is second order. We denote the 
two independent solutions by and 



An equally simple analysis for the d-dimensional case gives the following PDE:s: 

94 ^ dzkdzi dzk 2 

where k = 1, . . . ,d. In the d-dimensional case we have to distinguish between the cases 
when d is even or odd. The explanation for this is that when d is odd, the integrand 



of equation 1.5 has a branch point in the infinitely remote point, whereas in the even- 



dimensional case the eigenvalues Zi are the only branch points. To see this, set t = - and 



Zj = — . in the denominator of 1.5 



yj{t- zi)...{t- Zd) V(z - ai) . . . (z - ad) 
When z is close to the origin, w tends to 

w = re^""^'^, (2.6) 

which shows that the origin is a branch point in the z-plane if and only if d is odd. This in 
turn means that the infinitely remote point is a branch point for the integrand in equation 



1.5 if and only if d is odd. 

In order to get a single- valued integrand, we must cut the plane to prevent contours to 
encircle the branch points. This is done by pairwise connecting the branch points with 
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Figure 2: Branch cuts for the even-dimensional case. 



non-intersecting cuts. In the odd-dimensional case we also make a branch cut from one of 

the eigenvalues Zi to the point at infinity. This is illustrated in figure § and ^. 

In the even-dimensional case, any straight line not crossing any of the branch cuts is 



Z3 Zi Zs 
• • • 



Figure 3: Branch cuts for the odd-dimensional case. 



admissible if it is parallell to the real axis. In the odd-dimensional case, let be the 
branch point that is connected to the branch point at infinity. The admissble contours in 
this case are those that tend asymptotically to the lines Arg{zk) = and Arg{zj.) = vr and 
does not cross any of the other branch cuts. In this way we get two different classes of 
contours for the odd-dimensional case. The two classes are separated by the branch cut 
from Zk to the infinitely remote point. For the even-dimensional case there is only one class 
of admissble contours. The two independent solutions to equation 2A are denoted and 
J^, where is equal to zero in the even-dimensional case. 



3. Solution of the Ward identities 



3.1 Special cases 

Before we start with the general solution, we comment on some special cases. Consider the 
case when all the eigenvalues Zi are equal to one. This means that M2 = iMi. Moreover 
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the integral 



is just a constant so that 



_j2 12 

J,(l,...,l)= Idtjf-^ (3.1) 



Ad/2 ^ 

IdiSi, S2) = Jd oc , (3.2) 

^ ^ Vdet(Mi) V^et(Mi) ^ ^ 

which is a direct consequence of the following equation that was proved in [|l| 

^-S(xu...,x,)^^^ . . . OC / / {S{xi, Xd)) dxi... dxd- (3.3) 



Next, we consider the case when all the eigenvalues are equal, i.e. zi = . . . = Zd = z. In 



this case the differential equations 2.4 reduce to the Kummer-equation 



d'^Jd dJd d 

which has the solution 

Jd{z,...,z) = cU(-^,-,-ydM[-^,-,-y (3.5) 

where C/(a, 6, z) and M(a, 6, z) are the Kummer [/ and M functions respectively and c and 
d are constants. Notice that in one dimension this leads to the following well-known result 

where Ky{z) is the modified Bessel- function of the second kind. 
3.2 All eigenvalues distinct and non-zero 



In this section we solve the differential equations 2.4 when the eigenvalues of A are distinct 
and non-zero. The solution is found using Frobenius method, i.e. by assuming the function 
Jd can be expressed by an infinite series, 

00 

Jd= «n,.....4^+'^••4''"''^ (3.7) 

il,...,id=0 

where the ki are some real numbers to be determined. Since the integral defining the 
function Jd is completely symmetric in its arguments, it is clear that the coefficients ajj^,,,^i^ 
must be completely symmetric in their indices and that ki = k2 = ■ ■ ■ = kd, leading to 

CX) 

il,...,iti=0 
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Inserting this expression into equation 2A gives the fohowing algebraic equations 
_ i2 + k + l _ i3 + + 1 

Oil+2,i2,...,id — _^ _^ 2^11+1,12+1,13, ■■■,id _|_ _|_ 2'^*i+l'*2,«3+l,i4,...,id 

_ id + k + 1 _ ii + fc + 1/2 

i^ + k + 2''''+'''"-'''+' {ii + k + 2){h + k + if'^-'^' ^'^■^> 

and 

k{k -1) = 0. (3.10) 



Equation 3.10 is a consequence of the fact that any second order PDE will be a combination 
of two linearly independent functions. Note that in even dimension only k = yields a 
solution since the integral in this case is invariant under parity (i.e. the transformation 
Zi — > —Zi for all i). For k = we get, after some algebraic manipulations, the following 
expression for the coefficients 

(_l)(n+-+id)/2 (1) ... (i). 

an,...,^, = — r— TTZTTp ^^0' (^-^^^ 

il....id- U -h 2){n+...+ia)/2 
when ii + . . . ifi = 2n and zero otherwise. For k = 1 they become 

(_l)(n+...+..)/2(3)^^...(3) 



-Bo, (3.12) 



(n + l)!...(^. + l)!(f + 

when ii + . . . id = 2n and zero otherwise. In these formute, {z)n is the Pochhammer symbol 
defined by 

{z)n = z{z + l){z + 2)...{z + {n- 1)). (3.13) 
In summary we have the following expressions for the two independent solutions to equation 



2.4: 



and 



(-l)n(l) (I) 

4{zi, ■■■,Zd)=Ao Yl ■ ■ ■ 4^ (3-14) 

il+...+jd=2n *1 ^<i- U + 2Jn 



jO. ^ _ R det{M.. 

JdKZi, ■ ■ ■ ,Zd) - n\ 



det{Mi) 

where corresponds to /c = and J^^ to = 1. To simplify the notation we set = Q 
for d even. The integral I{Si, S2) is therefore given by 

h{Si,S2) = {J'd + 4) ■ (3.16) 

y'det{Mi) 

Taking the limit d ^ 00 only the terms corresponding to n = survive, leading us to 

Z^(5„5.) = ^=i=(A„ + ^fm). (3.17) 

^det{Mi) V det{Mi) J 
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3.3 All eigenvalues doubly degenerate and d even 

In this section we comment on the special case when d is even and all the eigenvalues of A 
are doubly degenerate. The differential equations in this case become 



,2 +Y.¥^ + -^-7r^ + JoD = ^, (3.18) 

ozt ^ ozkdzi dzk 



where the subscript DD stands for double degeneracy and /c = 1, . . . |. Notice that there 
are only | distinct eigenvalues in this case. We solve these equations in the same way as 
we did in the previous section, and the result is 

«l.--->«d/2=2n ^ ^ I 2^ I 

Let us look att the first few terms in the sum p. 19 . At the level n = there is just ^o- At 
level n = 1, we have the following 

2An 

^ {zl + ziZ2 + ziz^ + ...). (3.20) 



d + 2 

In general at level n, we have 



j ^dfn H2n{zi,...,Zd/2), (3.21) 



2 



where Hn{zi, . . . , z^/2) is the complete symmetric polynomial of degree n in the variables 
Zi. The expressions for the function Jj^f) can therefore be written 



/ 1\ " 1 
AD{zi,---,z,/2) = AoY,[--j JT^H'- 

V / n 



2n{zi,...,Zd/2), (3.22) 



leaving us with the following simple expression for the integral I {Si, S2) in this case 
3.4 S'0((i)-invariance 

Even though rewriting the integral Id{Si, S2) in terms of the eigenvalues of A = —iM2M^^ 
greatly simplifies the solution, it has the imediate drawback of hiding the 50(d)-invariance 
of Id{Si,S2)- To make this symmetry manifest we must replace the eigenvalues Zi with 
invariants of A = —iM2M^^. In principle this can be done via the Cayley-Hamilton 
theorem which (in one form) states that 

{t-zi)...{t- Zd) =t'^ + hf^-^ + ... + Id-it + Id, (3.24) 
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where Im are the invariants of A = —iM2M^ . Of course, this is an extremely difficult 
task except for the low-dimensional cases. In two dimensions we have 



Z2 



h + 
h- 



II + 4/2 



(3.25) 



so that 



i\j\n\ 



i,j=2n 



i! + 4/2 



(3.26) 



This expression gives the following manifestly 5*0 (2)-invariant formulation of /2('S'i) 'S'2): 

1 



4. Conclusion and Discussion 

In the present paper the non-Gaussian integral 



(3.27) 



Id{Si,S2) = / dxi 



(4.1) 



has been studied. A general formula for any dimension is given in equation 3.16. The 
formula is a function of the eigenvalues of the matrix A = —iM2M^^, where 



Si 
S2 



X Mix and 

X^M2X. 



(4.2) 



A formula for the infinite-dimensional case is given in equation 3.17, and depends only on 
the determinants of Mi and M2. This formula might be of use in quantum field theory. For 
two dimensions we give a formula with explicit S'0(2)-invariance, replacing the eigenvalues 
with the invariants of A. The formulas given in equations 3.14, |3.15| and 3.16 might be 
of help to calculate the integral discriminant J^|4, since any homogeneous action gives an 
inhomogeneous action when one of the coordinates is replaced by a constant, see equation 
(3) in [Q]. This has not been investigated in the present paper. 
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